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"ti ' Abstract 

The model under consideration is the two-dimensional Coulomb gas of ± charged 
hard disks with diameter a. For the case of pointlike charges (a = 0), the system is 

o ■ 

stable against collapse of positive-negative pairs of charges in the range of inverse 
temperatures < j3 < 2, where its full thermodynamics was obtained exactly [L. 
Samaj and I. Travenec, J. Stat. Phys. 101:713 (2000)]. In the present work, we 
£T) \ derive the leading correction to the exact thermodynamics of pointlike charges due 

to the presence of the hard core a (appearing in the dimensionless combination no , 
n is the particle density). This permits us to extend the treatment to the interval 
2 < (3 < 3 (the Kosterlitz-Thouless phase transition takes place at j3 = 4). The 
results, which are exact in the low-density limit no 1 — > 0, reproduce correctly the 
singularities of thermodynamic quantities at the collapse point (3 = 2 and agree very 
well with Monte-Carlo simulations. 
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1 Introduction and strategy 

The model under consideration is the two-dimensional Coulomb gas (2dCG), i.e., a neutral 

system of positive and negative unit charges qi = ±1 in a plane, interacting through the 

pair potential 

-qiqjInQri-Tjl/L), l^-r^Xr 

oo, |r» — r j | < cr 

Here, the logarithmic Coulomb potential is the solution of the 2d Poisson equation 
Af(r) = — 27r5(r). The Coulomb potential is regularized at short distance by a hard- 
core potential of diameter a around each charge. The system is studied as the classical 
one and in thermodynamic equilibrium, via the grand canonical ensemble characterized 
by the (dimensionless) inverse temperature (3 and the couple of equal particle fugacities 
z + = z_ = z. The fugacity z has dimension [length] -2 . Within the grand canonical 
formalism, the length scale L in (|1.1| ) manifests itself as the rescaling of z, z — > L^l 2 z. We 
shall set L to unity for simplicity, keeping in mind that the true dimension of the rescaled 
z is [length]' 3 / 2-2 . The corresponding particle number densities n + = ra_ = n/2 (n is the 
total particle density) enter into the formalism in the dimensionless combination na 2 . 

For small values of the dimensionless density na 2 , the famous Kosterlitz-Thouless 
(KT) transition PJ of infinite order takes place at a specific density-dependent inverse 
temperature /3kt 0, ||- I* 1 the high-temperature conducting phase (3 < /3kt> the ef- 
fective potential between infinitesimal external charges decays exponentially due to the 
perfect screening by the positive and negative charges of the Coulomb system. In the 
low-temperature dielectric phase (3 > /3kt, the system charges form dipoles and no longer 
screen an external charge, so that the effective potential between infinitesimal external 
charges is proportional to the bare logarithmic potential. At high enough na 2 , the KT 
critical line splits into a first order liquid-gas coexistence curve (for Monte-Carlo (MC) 
simulations, see refs. 0, 0], for theoretical computation, see ref. 0). 

In the low-density limit na 2 — > 0, which is of special interest in general and also in this 
paper, the thermodynamic behavior of the 2dCG as a function of /3 undergoes fundamental 



changes at two points: j3 c = 2 (the collapse of pointlike particles) and /3kt = 4 (the KT 
phase transition). The first - collapse - point reflects the fact that, for the case of strictly 
pointlike particles a = 0, the singularity of the Coulomb potential v(r) ( |1 . 1|) at the 
origin prevents the thermodynamic stability against collapse of positive-negative pairs of 
charges (or, equivalently, the corresponding Boltzmann factor r _/3 is not integrable at 
short distances in 2d) for (3 > 2. Thus, for < j3 < 2, the system of pointlike particles 
is thermodynamically stable and the introduction of a hard core around particles is a 
marginal perturbation which does not change the thermodynamics substantially. On the 
other hand, for 2 < /3 < 4, the introduction of a hard core is inevitable for avoiding 
the collapse: when one calculates thermodynamic quantities and at the end takes the 
limit a — > (with z being fixed), while the density, the free energy and the internal 
energy per particle diverge due to the collapse phenomenon, the specific heat and the 
truncated (Ursell) correlation functions are expected to remain finite ||. In spite of the 
tendency to the collapse into neutral pairs of charges, there still exist free charges which 
are able to screen and the system remains in its conducting phase up to the KT phase 
transition at point /3kt = 4. We would like to stress that, for a given fugacity z and when 
(3 > 2, although n — > oo as a — > 0, the dimensionless density goes to the limit of interest, 
no 2 -»■ 0. 

In what follows, we shall summarize in detail the known results in the two qualitatively 
different regimes (0 < (3 < 2) and (2 < (3 < 4), and at the collapse point /3 C = 2. 

In the stability range of inverse temperatures < /3 < 2, as has been already men- 
tioned, the thermodynamics of the 2dCG is well defined even for the case of pointlike 
particles, a = 0. The density derivatives of the Helmholtz free energy, like the pressure p, 
can be calculated exactly by using a simple scaling argument. For instance, the equation 
of state 

(3p = n(l-^ (1.2) 

has been known for a very long time 0. The temperature derivatives of the Helmholtz 
free energy, like the internal energy U or the constant volume (surface in 2d) specific heat 



Cy, are nontrivial quantities, the calculation of which can be based on an explicit density- 
fugacity relationship. This relationship was obtained only recently || via a mapping of 
the 2dCG onto a classical 2d sine-Gordon theory with a specific normalization of the 
cos-field, and then by using quite recent results about that integrable field theory || [l(J . 
Explicitly, 
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where T stands for the Gamma function The density-fugacity relationship Q1.3Q was 

checked on a few lower orders of its high-temperature /3-expansion by using a renor- 

malized Mayer expansion in density, valid just in the stability regime. For fixed z, the 

particle density exhibits the expected collapse singularity as (3 — * 2~: 

4nz 2 , N 

n~—p (1.4) 

This behavior can be derived by using an independent-pair picture of the system around 
the collapse point ||, which is another check of the exact results. Based on the density- 
fugacity relationship (|1.3| ), the complete thermodynamics of the pointlike CG can be 
obtained by elementary means in the whole stability interval < (3 < 2 [H. 



At the collapse point f3 c = 2, by the continualization of Gaudin's lattice model fll 
which is expected to have the same properties as the 2dCG in the low-density limit, the 
truncated many-body densities (Ursell functions) were found in refs. |12| , p~3|j . These 
densities have the remarkable property of going to well-defined limits as na 2 vanishes (as 
is believed, this property lasts up to the KT phase transition), identical to the densities 
of an equivalent Thirring model at the free-fermion point. The knowledge of all truncated 
many-body densities at /3 C = 2 for na 2 = permits one to extract the leading parts of 
thermodynamic quantities at /3 = 2 which, for a fixed fugacity z, do not vanish in the 
low-density limit na 2 — > |T2"||. Namely, 



n = Attz 2 
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KairzJ J 8 ' \ln(cnrz) 
Here, u cx = (E)/N is the excess (over ideal) internal energy per particle, dy = Cy^/N is 
the excess specific heat at constant volume per particle, and C is the Euler constant. 

The region of inverse temperatures 2 < j3 < 4 is usually studied by using the Mayer 
series expansion of the specific grand potential in fugacity. It was proven that each term 
of the ^-series converges in the insulator region (3 > 4 [[14], [15]]. For (5 < 4, the existence 
of infinitely many thresholds at inverse temperatures 
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lying between j3\ = (3 C = 2 and /?oo = /3kt = 4, was observed: if j3 > fy, only the Mayer 
series' coefficients (cluster integrals) up to the order 21 are finite, and the cluster integrals 
of order > 2/ exhibit the large- distance divergence in the infinite- volume limit. Below the 
collapse point (3 < 2 (where the density format is appropriate) all cluster integrals diverge. 
The free energy is supposed to have at points {j3i}fZi a logarithmic dependence on the 
cut-off (in our case na 2 ) [|T^]. Points {/3i}^ 1 were conjectured to correspond to a sequence 
of transitions from the pure multipole insulating phase ((3 > 4) to the conducting phase 
(/3 < 2) via an infinite number of intermediate phases ||14|| . 

Such a conjecture was later denied by Fisher et al. [I7|. In the limit na 2 — > 0, they 



proposed the following ansatz for the equation of state in the fugacity format :[] 
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where the temperature-dependent exponent is 
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1 we rescale z by cr' 4 @'> 2 in the ansatz in order to work with dimensionless quantities 



The function e vanishes as na 2 — > 0. At points A [see relation ( |1.6|) 1, V(A) = I- The 
coefficient b^((3) and the function e^o^ 4- ' 3 )/ 2 , (3) were suggested to be analytic in the whole 
conducting regime < f3 < 4. The coefficients {b2i(P)}f^i were suggested to be zero for 
< /3 < 2, and finite, analytic in /3, for 2 < /3 < 4. The divergence of the coefficients of 
an analytic expansion in z 2 was related to the appearance of the anomalous term b^z 2 ^ in 
(|1.7[). Moreover, at Act = 4, iP(/3kt) becomes infinite and the anomalous term disappears 
from the ansatz, in agreement with the results of refs. |T^, I3| . The findings of this work 



support the ansatz (|I~T|), however, the supposed analytic behavior of the coefficients b^((3) 
and {b2i{(3)}^Zi in the whole interval of (3 G (0,4) is apparently wrong. For example, in 
the stability region < (3 < 2 and in the limit of pointlike particles na 2 = 0, the equation 
of state ( |OD , when combined with the ansatz (|1.7|) with e = 0, yields 



W) = (i - f ) ^bi ( L9 ) 

Substituting n/z 4 ^ 4 '^ by the rhs of (|T3|), we get the exact b^((3) oc tg(vr/?/[2(4 - /?)]). 
This function has the simple pole just at the collapse point (3i = (3 C = 2, and its analytic 
continuation to the region 2 < (3 < 4 has the simple poles at {j3i}fZ 2 given by ( |1.6|) (see 
also the next paragraph). 

The aim of this paper is to derive the leading correction to the exact thermodynamics 
of pointlike charges due to the presence of the hard core of diameter a, and in this way to 
extend the treatment of the 2dCG beyond the collapse point f3 c = 2. First of all, we have 
to understand the analytic structure of a "naive" continuation of the exact results for the 
particle density of pointlike charges ( |1.3| ) to the range of inverse temperatures 2 < (3 < 4, 
in order to see which singularities should be removed by the consideration of a hard core. 
The only source of singularities in n ( |1.3| ) is the tg-function whose decomposition into 
simple fractions [0 gives 



tg(7r/?/[2(4 -/?)]) 2^, 1 (4-/3) 2 (4-A) 2 ( , 

tt/3/[2(4 - /?)] " ^ £j (A - (3) (3(4 - A) + A(4 - P) l ' } 

The thresholds {A}^i> given by ( |1.6[ ), manifest themselves as simple poles at which 

n exhibits the discontinuity from linig^- n((3) — > oo to linig ~+ n(j3) —>■ — oo. The 
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introduction of the hard core a must remove the artificial singularities of n at {/3i}fZ\ and 
must make n positive as is required by its basic definition. 

We reach our aim by combining the zeroth-moment (electroneutrality) sum rule for 



the pair charge-charge correlation function [19], valid in the whole fluid regime < (3 < 4, 



with the short-distance expansion of that charge-charge correlation. The short- distance 
expansion of the charge-charge correlation can be done, in principle, systematically. Every 
new term of order "Z" should remove the singularities of the pointlike n at j3i and makes 
the particle density positive up to /?j+i. We were able to establish the leading term of 
order "1" which permits us to extend the treatment of the 2dCG to the interval 2 < /3 < 3. 
The advantage of our approach with respect to the standard ones based on the fugacity 
expansion is that, in the whole interval < (3 < 3, our estimates are exact in the low- 
density limit na 2 —>■ 0. Moreover, the leading parts of thermodynamic quantities at j3 = 2 
( |1 . 7|) are reproduced correctly. The agreement with the MC simulations @] is very good. 
The paper is organized as follows. The method is described in Section 2. The complete 
thermodynamics is derived and then tested in various limits in Section 3. In section 4, the 
comparison of our results is made with the MC simulations and some concluding remarks 
are given. 

2 Method 

In the conducting phase of the two-component plasma with the Coulomb plus an arbitrary 
short-range pair interaction of particles, the zeroth and second moments of the charge- 
charge density are determined exclusively by the long-range tail of the Coulomb potential 



19]. The starting point of our calculation is the zeroth-moment (electroneutrality) sum 
rule 

^ = /[^,-,(r)-^(r)]d 2 r (2.1) 

Here, with the notation n q (r) = Y^t^q^ij — r «) f° r the microscopic density of particles 
of charge q — ±1 at position r, n g = (n q (r)) = n/2 and the Ursell function U g>q i(r, r') = 



,(2)T 



Uq } q>(\r — r'|) (denoted as p , in refs. Jl2|, [13[]) is denned by 



C/„„/(r,r') = (n,(r)7v(r / )) - n,5g ig /5(r - r') - n g n g / (2.2) 



For the Coulomb system of interest (|1.1|), the difference U q - q — U q>q vanishes inside the 
hard core, and the total particle number density n is given by 



oc 



n(z, a) = 2 / 2wrdr [U q _ g (r; *, cr) - f7 M (r; z, a)] (2.3) 

Hereinafter, we omit in the notation the dependence of quantities on f3. 

The dependence of the density n on a in Q2.3Q comes from the cutoff in the integration 
over r and from the a-dependence of the Ursell functions themselves. These Ursell func- 
tions are supposed to be well defined and finite for the zero density, no 2 = 0, in the whole 
conducting regime < (3 < 4, including the collapse interval 2 < j3 < 4. This belief is 
strongly supported by the finite values of the Ursell functions at the collapse point /3 C — 2 



12, m-. 



' 2\ 2 



U q ,- q (r;z,0) = 1 — 1 ^f(mr) (2.4a) 



' 2\ 2 



l/ M (r;z,0) = - I — J K^(mr) (2.4b) 

with fixed m = 2ttz, Kq and K\ are modified Bessel functions. We can thus write 

U qtq/ (r;z,cr) = U m >(r;z,0) + A qjql (r;z,cr), r > a (2.5) 

which defines A q>q /(r; z, a), vanishing when a — > 0, as the change of U q>q /(r; z, 0) due to 
the introduction of the hard core a < r to pointlike particles. Subtracting Eq. ( [2.3|) with 
a > and the same equation with a = 0, one arrives at 

n(z, a) - n{z, 0) = -2 f 2vrrdr [E/ tf ,_,(r; ^, 0) - U q , q (r; z, 0)] 

Jo 

("OO 

+2/ 27rrdr[A gi _, 2 (r;2;,a)-A gi(? (r;2;,(T)] (2.6) 

Since n(z, 0) is defined only in the stability regime < (3 < 2, for the time being we shall 
restrict ourselves to this range of (3. We now make a heuristic assumption analogous to 
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that made at the point (3 C = 2 in refs. [|T|, [H|: in the low-density limit na 2 — > and for 
< (3 < 4, one can neglect the quantities A q ^± q in Eq. ( j2.6|) . Consequently, 

n(z, a) = n{z, 0) - 2 f 2vrrdr [J7,,_,(r; z, 0) - U m (r; z, 0)] (2.7) 



Our assumption is equivalent to saying that in equation (|2.6| ) only the contribution oc 
Jq rdr[U q _ q (r; z,0) — U q>q (r; z,0)] with a > is enough for removing, via a systematic 
short-distance expansions of U q ^± q (r; z,0), term by term the singularities of n(z, 0) at 
points {(3i}fZi (see the Introduction). This scenario will be verified at the first singular 
point Pi = j3 c = 2. 

The next step is to construct the short- distance expansions of U qy ± q (r; z, 0) in the 
integral on the rhs of ( |2.7[ ). For small enough (3, the short-distance expansion of the 
Ursell functions is dominated by the Boltzmann factor of the corresponding pair Coulomb 



potential 0, £| • I* 1 the case of oppositely charged particles, one has 

U q _ q (r; z, 0) ~ z?r~ as r ->• (2.8) 

valid in the whole interval < (3 < 4. Note that since Ki(mr) ~ l/(mr) for r — » 0, 
£/g_ g (r; 2,0) at /3 C = 2 [see relation (|2.4a|) 1 satisfies (j2.8|). In the case of the same charges, 



the leading term has a more complicated structure, 



r? for < 3 < 1 
C/,, g (r;2,0) oc <( as r -> (2.9) 

r 2 ~P for 1 < 8 < 2 



The change of the power-law behavior is caused by the divergence of the prefactor to r at 



(3=1 [pi]] . Considering in ( p.4b| ) that K (mr) ~ — lnr for r — > 0, the logarithmic behavior 
of U q:q (r; z, 0) at Q c = 2 can be understood as a limiting case of r 2_/3 in ( p.9|) . Within the 
sine-Gordon representation of the 2d pointlike CG, the formula corresponding to ( |2.8| ) 
is known as the conformal normalization of the cos-field. For such a theory, the short- 
distance expansion of correlation functions is available by using the Operator Product 



Expansion [521 , as was explicitly done in ref. [23]. Although this method allows one to 



construct systematically the short-distance expansions of U q: ± q (r; z, 0) with coefficients 

9 



expressed in terms of Dotsenko-Fateev integrals p4] , it applies only to small values of f3 



and does not describe, e.g., the change of the behavior at (3 = 1 (2.9). In any case, in the 
interval < j3 < 4, the charge-charge combination of the Ursell functions is dominated at 
short distances by ( |2.8|) , 

E/, _ g (r; z, 0) - U M (r; z, 0) ~ z 2 r^ as r -► (2. 10) 

and we shall consider just this leading term. 

Inserting ( |2.10| ) into ( |2.7p , one obtains the basic formula 



a 2 ' 13 
n(z, a) = n{z, 0) - 4vrz 2 — - (2.11) 

Although this result was derived in the region < (3 < 2 where n(z, 0) is well defined 
by ( |1.3|) , it is reasonable to assume the continuation of ([2.1 1| ) beyond the collapse point 



j3 c = 2: both the sum rule ( |2.3| ) and the leading short- distance expansion ( [2.10|) , which 
play the crucial role in the derivation of ( |2.11| ), remain valid up to the KT transition 
at /?kt = 4. As will be shown in the subsequent section, the leading correction term in 
Eq. ( [2. 11|) removes the singularity of n(z, 0) at (3\ = (3 C = 2, and provides an adequate 
description of the 2dCG with na 2 small up to fa = 3 where another singularity of n(z, 0) 
occurs. This singularity at fl% — 3 should be removed by the next term of the short- 
distance expansion ( |2.10| ) inserted into ( |2.7| ), however, as was mentioned above, it is not 



simple to get the explicit form of that term for such a large value of (3. 



3 Thermodynamics 



Based on the density-fugacity relationship ( [2.11| ), we derive in this section the thermody- 



namics of the 2dCG in the restricted region of interest < j3 < 3. Our findings will be 
compared with the known results and conjectures reviewed in the Introduction, namely 



with the exact formulae at (3 C = 2 (|1 . 5|) , with the predictions based on the independent- 



pair collapse picture for f3 > 2 and in the limit o — > |J, and with the proposal (|1.7| ) for 
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the equation of state. The stability regime < (3 < 2 will not be discussed in detail since 
there the introduction of the hard core to pointlike particles is a marginal perturbation. 



We first express relation (|1.3|) for the density of pointlike particles n(z, 0) in a more 
convenient form, 



,,(,« = ??*» 



where we have introduced the function 

/J/(4-/3) r 
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(3.2) 



(1/tt)P (1/2 + #[2(4 - /?)]) vr/?/[2(4 - (3)\ 

With regard to the collapse singularity ( |1.4| ), $ was choosen such that $(/? = 2) = 1. The 
Taylor expansion of In $(/3) around (3 = 2 thus reads 



ln$(/?) 



(hi7r + C)(/3-2)H-~(hi7rH-C)(/3 



+- 



lnvr + C-^C(3) 



(/3-2) 3 + -- 



(3.3) 



The function 3>(/3) is positive in the interval < (3 < 8/3, it crosses zero at point (3 = 8/3 
(which is not exceptional from any point of view) and diverges to — oo as j3 —>■ 3. Using 
the representation (|3TT| ) in Eq. (|2.11| ), one gets for n = n(z, a) 



n 



An 

2^p 



4/(4-/3) 
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(3.4) 



Notice that this relation is dimensionally correct - it is expressible in terms of dimension- 
less quantities 

f = az 2 ^~V (3.5) 



and na 2 , or the packing fraction 



77 = -una 
4 



(3.6) 



used in the MC simulations 0. Taking into account the expansion (|3.3|) in Eq. ( |3.4|) one 
sees that, for fixed z and nonzero a, the density n is finite at (3 C = 2; its value coincides 
with the expected result ( |1.5a| ). Formula (|3.4| ) can be analytically continued to the region 
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2 < (3 < 3. For (3 > 2 and small a, the term [crz 2 ^ 4- ^] 2- ^ becomes larger than $(/?) and, 
when combined with the negative denominator (2 — (3), it implies the positive sign of n 
up to (3 = 3. When a — > 0, this term and consequently n diverge in the region 2 < (3 < 3 
as a consequence of the collapse phenomenon. 

We proceed by the derivation of the equation of state. The grand canonical potential 
Q is determined by the thermodynamic relation 

n=z S ±pn (3.7) 

oz 

where V(— ► oo) is the volume (in 2d, the surface) of the homogeneous system. With 
respect to the boundary condition (—(3VL/V)\ z= q = 0, the integration of (|3.7|) , with n 
substituted from Eq. (|3.4|), results in the equation of state for (3p = —j3Q/V, 

(3p = ^f|^(/3)^(^) - ^|^y (za^f (3.8) 

Eq. (|3.8|) has the form of the ansatz ( |1.7|) with 

h ^ m , 5, = -^ ,3,) 



In contrast to the conjecture made in ref. [17|, these coefficients are not analytic functions 



of (3 at (3 — 2. Let us now rewrite the equation of state ( |3.8|) with the aid of the density- 
fugacity relationship fl3T4|) as follows 

^- l | ( 2 -^ (3 10) 

where ^ is defined by ( |3.5|) . At (3 C = 2 and for nonzero a, (|3.10|) agrees with relations 
( |1.5a| ) and (|1 . 5b| ) . In the limit of pointlike particles a — » 0, our equation of state (|3.10|) 



reduces to the one obtained by Hauge and Hemmer 0: in the stability regime < (3 < 2, 
the hard-core correction £ 2-/3 is negligible for small a and the a — > limit corresponds to 
(|1.2|), while for 2 < /3 < 3 £ 2 ~^ diverges when a — >• and so j3p/n — > 1/2, i.e., the pairs 
of ±1 charged particles collapse into neutral "free" particles of half density n/2. 
Finally, the (excess) dimensionless specific free energy / = (3F e */N reads 

ft R\ P Q ^ 

f{n,(3) = — - + \nz 
nv 

= -I- f-^+ln, (3.11) 

2 4($-^ 2 -' 3 ) v ; 
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where the implicit dependence of the fugacity z(n, j3) on the particle density n is deter- 
mined by Eq. (|3.4j) . According to the elementary thermodynamics, the (excess) internal 
energy per particle u ex is given by w ex = df(n,fl)/d(3, explicitly 



ex 1, 1 1 

u = - In a 



2-|W + $ln£ 



2 2(2 - (3) 2($ - £ 2 -/3 

and the (excess) specific heat at constant volume per particle dy is given by dy/kg 

-(3 2 d 2 f{n,(3)/d(3 2 , explicitly 

J_£E 1 (4 - /?)$" 

~]PY B ~ 2(2-/3) 2 + 4($-e~ /3 ) 

1 



(3.12) 



2(2 -/?)($- f 2 -0) [$ - (2 - /3/2)£ 2 H*] 
x J (2 -/?)$' 



$+h2-|W 



($-^)l 



+ (2-/3)^(lnO[(l + lnO* + (4-/3)$'] (3.13) 

At /3 C = 2 and for nonzero <r, w ex is identical to ( |1 . 5c|) , while 

c e v (g = 2) 1 [ln(Tr^) + C] 3 + 2[ln(7r^) + Cf 17 1 

fc B 3 l + 2[ln(vr^)+C] 12 U ' [\n{-Kza) + C] l J 

The expansion of this expression into the Laurent series in l/[ln(7T2(j) + C] reproduces 

the leading terms in ( |1.5d| ). In the region 2 < (3 < 3 and in the limit a — > 0, u ex diverges 

to — oo due to the term (lna)/2 involving the energy of the collapsed pair of particles, d^ 

behaves differently: £ 2-/3 becomes large for small a and kills all terms on the rhs of Q3.13 ), 



except of the first one. The consequent finite result Cy/ks = /3 2 /[2(2 — /3) 2 ] coincides 
with the finding of ref. ||. 

4 Comparison with MC simulations and conclusion 

In the previous two sections, we have derived the thermodynamics of the 2dCG by adding 
the leading hard-core correction term to the exact density-fugacity relationship for the 
Coulomb system of pointlike particles. We have successfully tested various limits of our 
results. In this section, the comparison is made with the MC simulations of the 2dCG Q. 
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Figure 1: The plot of Cy/ks vs. the packing fraction r\ at (3 C = 2: present result (|3.14j) 
(solid line), formula (|1.5d| ) (triangles) and MC simulations (circles). 

The results for the heat capacity dy at the collapse point f3 c = 2, as the function of 
the packing fraction rj (|3.6| ), are presented in Fig. 1. Our formula ( |3.14|) , represented by 
the solid line, provides the estimates of Cy* consistent with the MC data (circles) within 
the error bars. The agreement with the MC simulations is much better than in the case 
of the previously derived formula (|1.5d|) (triangles), which involves the first three terms 
of the Laurent expansion of ( |3.14| ). 

In Figs. 2 and 3, we represent by solid lines our plots of the internal energy u cx ( |3.12| ) 
and the heat capacity Cy ( |3.13j ), respectively, versus the inverse temperature (3, for two 
values of the packing fraction rj = 5 x 10~ 4 and 5 x 10~ 3 . The agreement with the MC 
data (circles) is very good, but getting worse when approaching point j3 = 3: here, the 
next (neglected) term of the hard-core corrections, which eliminates the next singularity 
of n(z,0) at /?2 = 3, starts to be important. 

In conclusion, the extension of the thermodynamic treatment of the 2dCG around 
and beyond (3 2 = 3, up to /3kt = 4, requires to construct a systematic short- distance 
expansion of the Ursell functions for pointlike charges. It is clear that the relevant hard- 
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Figure 2: The plot of w ex vs. the inverse temperature (5 for 77 = 5 x 10 4 and 5 x 10 3 : 
present result ( |3.12|) (solid lines) and MC simulations (circles). 
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Figure 3: The plot of Cy/hs vs. the inverse temperature (3 for rj = 5 x 10 4 and 5 x 10 3 : 
present result ( ^.13| ) (solid lines) and MC simulations (circles). 
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core corrections must subtract all simple poles of the pointlike-particle density n(z, 0) 

( |1.3|) at {Pi}fZ\ given by ( |1.6|) . Thus, the density- fugacity relationship should be of the 

form 

n{z,*)=n{z,Q)-£f@ i (*>o*-0) 1 (4.1) 

1=1 Pi ~ P v 

where the coefficient c 2 i(P) (I = 1,2, . . .) is regular at (3 = fa. According to our result 
( [2. 11| ), C2{(3) = Ait. Taking into account the formula (|1.3|) for n(z, 0), the exponent of z, 
4/(4 — j3), is equal to 21 at the pole /3 — fa and the exponent of o is put by dimensional 
reasons. Note that since n = zd(/3p)/dz, (|4.1| ) automatically leads to the ansatz (|1.7|) 
with e = [because we have not taken into account in ( |2.6|) "slight" hard-core changes 
of the correlation functions Aq j±(? defined by fl2.5Q]. The coefficients b^ and {&2«}^i are 
singular functions of fa The Ith term in the sum on the rhs of ( f4.1| ) has the correct 
behavior in the limit a — > 0: it goes to when (3 < fa and removes the singularity of 
n(z, 0) at P = fa (and simultaneously gives rise to a logarithmic dependence on the hard 
core at this point). A systematic generation of the coefficients q(/3) is our task for the 
future. 
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